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X. 

^ . Abstract 



We study the energy dependence of collective (hydrodynamic-like) nuclear 
matter flow in 400-1970 A MeV Ni+Au and 1000-1970 A MeV Ni-hCu reac- 
tions. The flow increases with energy, reaches a maximum, and then gradually 
decreases at higher energies. A way of comparing the energy dependence of 
flow values for different projectile-target mass combinations is introduced, 
which demonstrates a common scaling behaviour among flow values from dif- 
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fererent systems. 

PACS numbers: 25.70.Pq, 25.75.Ld 



Typeset using REVT^ 



The study of nuclear matter over a wide range of temperatures and densities and the 
determination of the Equation of State (EOS) of nuclear matter continue to be of consid- 
erable interest Lacking the possibility of comprehensive studies of nuclear matter in 
bulk (as in neutron stars), one resorts to the study of transient, finite systems provided by 
nucleus- nucleus collisions (over a wide range of energy). It is now clear that the extraction 
of EOS information from nuclear collisions requires a comprehensive set of measurements of 
collision observables, which can be compared to realistic microscopic calculations involving 
the nuclear matter variables. 

At GeV per nucleon energies, where the collision velocity exceeds that of nuclear sound, 
the collisions produce densities several times higher than ground state densities and exhibit 
compression-induced flow of nuclear matter [^. However, it was not until the analysis of 
events from the Air Plastic Ball/Wall 0] and Streamer Chamber detector systems at the 
Bevalac that the matter flow characteristics could be studied and quantified for a range of 
systems and energies. More recently, at the Bevalac at LBNL, the EOS Collaboration carried 
out a comprehensive set of measurements over a wide range of energies, and projectile- 
target combinations [^. Similar studies are underway at GSI in Darmstadt and, with lighter 
projectiles, by the DIOGENE Collaboration at Saclay EOS Collaboration data has been 
used to study flow for the Au + Au system at lab energies ranging from 250 to 1150A MeV 
. Particle flow for protons, deuterons and alpha particles has been determined, using the 
transverse momentum method 0. The flow is found to increase with particle mass A, and 
with energy up to projectile energies of 1150A MeV where it tends to level out at values 
close to the Plastic Ball data 0. 

Here we present analyses of recent EOS Ni + Au data with energies between 400 and 
1970A MeV and EOS Ni + Cu data at energies of 1000, 1500 and 1970A MeV. This study 
with the Ni beam allows the use of higher energy per nucleon projectiles, and so extends the 
energy of our flow measurements beyond the 1150A MeV limit of the EOS Au + Au data 
10. We also present a comparison of the flow values with predictions of a BUU model from 
Bauer et al. 0, and introduce a scaled flow which allows the comparison of flow data from a 
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variety of projectile-target mass systems. In total, the data (supported by the calculations) 
provide the strongest evidence yet that, with increasing energy, flow reaches a maximum 
near 1000 A MeV, and then declines. 

The EOS Collaboration detector systems have been described in Refs. 4 and 6. The 
data presented here were obtained using the EOS Time Projection Chamber [|T^], situated 
in the magnetic fleld of the HISS Magnet. The TPC provides fairly unambiguous particle 
identification, as well as a measurement of momentum, for particles of charge up to Z = 
6. Particle ID is ambiguous for rigidities above 2.4 GeV/c. Thus, some misidentification 
will occur and its effect is represented in the uncertainties. The target is just upstream 
from the TPC and this results in a large and nearly seamless acceptance. Laser beam 
calibrations provide a check of the corrections for B field inhomogeneities, and measurements 
of drift velocity. Simulations have been performed to study the geometrical acceptance of 
the detector and to provide acceptance corrections. 

In the present experiment, the thickness of the Au target was 730 mg/cm^ correspond- 
ing to about 1% interaction probability. The trigger was provided by a scintillator just 
downstream from the target. The MUSIC (Multiple Sampling Ionization Chamber) detec- 
tor, downstream from the TPC, provided an on-line check on the threshold charge of the 
heaviest fragment allowed by the trigger. For this analysis, data for a wide range of impact 
parameters were taken. We used charged particle multiplicity as a measure of the collision 
centrality, and have adopted the Plastic Ball @| convention by dividing the events into five 
multiplicity bins with bin Ml corresponding to the most peripheral and bin M5 having the 
most central events. For the Ni + Au 400, 600, 1000 and 1970A MeV the numbers of events 
analyzed were 32k, 26k, 33k and 26k respectively. The Ni + Cu 1000, 1500 and 1970 A 
MeV analyses were done on recently analyzed data sets of respectively, 35k, 19k and 48k 
events. For each event, the reaction plane was determined using the transverse momenta 
of the particles, as proposed by Danielewicz and Odyniec [0. The plane is determined by 
the vector Q = J2i WiP \ and the incident beam direction. Here p \ is the transverse mo- 
mentum of particle and Wi is a weighting factor defined to maximize the contribution of 
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high rapidity particles to the Q vector determination. Normally, for symmetric collisions of 
equal A, Wi is taken to increase linearly with rapidity up to Wi = ±1 at ly'c^l = S and to be 
constant at ±1 for \y'c^\ > S. Using the sub-event method of Ref. 7 we found the optimal 
reaction plane determination to be made with S ~ 0.7. We also found that varying 5 by 
0.1 results in only a 1% change in the flow values. For the asymmetric ^^Ni + ^^"^Au, the 
nucleon-nucleon and nucleus-nucleus center of mass frames are different. There is therefore 
some ambiguity as to which reference frame to use for determining the weighting factors. 
We went through an iterative procedure as follows: We start out in the nucleus-nucleus cm 
frame, do the Q weighting in that frame, and then project the pt of each particle onto the 
reaction plane. We refer to this projection as for each particle in the event. Plotting 
the (px/A) (henceforth referred to as Px) vs. rapidity for all events at a given energy then 
yields the typical S-shaped curves. These S curves do not cross the Px = axis at i/c — 
where i/c here refers to the rapidity in the chosen frame, in this case the nucleus-nucleus cm 
frame. We next look at Px vs. lab rapidity and note where the curve crosses the Px = 
axis. This crossing rapidity defines the velocity of the new frame in which we do the next 
weighting. A couple of iterations yields a stable crossing value of y. Finally, we need to 
correct px for the fact that we project onto an imperfectly known reaction plane. For this 
we use the sub-event method of Ref. 7. These corrections increase the Px values from 10% 
(for the more central Ni + Au events) to 20% (for the more peripheral Ni + Cu events). 

Figure 1 shows the S-shaped plots of Px vs. y' for the four Ni -\- Au systems and for 
multiplicity bin M4. Here, for the moment, we include only protons, for which we have 
unambiguous particle identification. Following the Plastic Ball analysis we define the flow, 
F, as the slope {dpx/dy')p^=o at the zero crossing (generally around y' — 0.35). The slope is 
calculated from a linear fit to the data. Similar F values are obtained using a cubic fit. The 
imperfect asymmetry of the S-curve with respect to the Px — ^ axis is due to both the lower 
acceptance in the lower rapidities, as well as reflecting transverse momentum conservation 
in the asymmetric collision: The larger number of (mainly target) particles (both spectator 
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and participant) at low rapidity, is balanced at higher rapidity, by fewer particles having 
larger average pt values. From the geometrical acceptance studies we found that acceptance 
corrections on the flow values are on the order of 5% for the lower energy Ni + Au 400A 
and 600A MeV systems . Flow value corrections above 600A MeV are negligible since the 
detector has full acceptance for the region where the fitting is performed. 

In Figure 2 the extracted values of the slope, F (for protons), at y'{px = 0) are plotted 
vs. projectile kinetic energy per nucleon and compared to the flow predicted by a basic BUU 
model 1^. This model does not include composite particle formation so the comparison is 
with the experimental proton F values. With the inclusion of bound protons (from d, t, 
etc.) the trend of the F values, as a function of energy, is similar. It should be noted 
that the BUU model used here does not include momentum dependent interactions. It is 
expected that adding these will increase the predicted BUU flow values, bringing them to 
better agreement with the data. The hard EOS BUU flow values appear to reach their 
maximum at a lower energy than the soft EOS. In this aspect the hard EOS BUU is closer 
to the energy dependence of the data. However, neither case gives quantitatively the energy 
dependence, and work on this aspect continues [O]. The main point here is that the BUU 



model does predict a rise and fall of flow. The energy dependence in Fig. 2 is consistent 
with the EOS Au + Au proton flow data of Partlan et al. ^ and the Plastic Ball data p| 
which show the flow beginning to plateau above 800A MeV. 

We have studied other flow observables: the maximum value of pr^, the maximum value 
of F, regardless of the corresponding y value, and the sum J2Px, summed over all values 
for y' > y'{px = 0). All of these observables give a flow energy dependence consistent with 
that for F. We have also analyzed the data using the reaction plane independent flow signal 
quantity ||I2[ proposed by the FOPI group. In addition we have examined the flow angle 
dependence on multiplicity for the 4 energies. Both the flow angle and the FOPI signal show 
an energy dependence (rise and fall) similar to that of F in Fig. 2. 

It can be argued that flow should be determined from the laboratory rapidity, y, rather 
than from reduced rapidity, y' = y/yhcam- Then flow becomes Fy = dp^/dy = F/yp where 



Up is the projectile beam rapidity in the lab frame. Plotting Fy vs. energy we find that the 
rise of fiow is reduced while the decline becomes more significant. 

The rise and fall of fiow is explicable or at least qualitatively reasonable. As energy 
increases, the nucleon-nucleon cross sections become more forward peaked. The mean 
transverse momentum at first rises rapidly with energy and then is almost constant above 
P2 — 2 GeV/c. Thus, fiow as a manifestation of multiple scattering effects should eventually 
fall as pz continues to increase. 

It is of great interest to compare the flow values for a wide range of data. To allow for 
different projectile-target {Ai,A2) mass combinations, we divide the flow value by {A\^^ + 
A^2^) and call Fg = F/{Ay^ + A2^^) the scaled flow. In recent calculations of flow ||13|, the 
authors use an F/A^^^ scaling for symmetric systems, and argue that, for a given energy, 
the flow should scale with collision (compression) time. This suggested the {A\^^ + ^2^'^)^^ 
scahng used here. Figure 3 shows a plot of Fs vs. energy per nucleon of the projectile. We 
include here data from the EOS experiment (solid points), along with values derived from 
other experiments [P|,P^|-P^ for a variety of energies and mass combinations. As closely as 



possible the data selected correspond to Plastic Ball multiplicity bins M3 and M4 or to an 
equivalent range of impact parameters. For the EOS and Plastic Ball data all the isotopes of 
Z = 1 and 2 are included, except for [|l6l where the data is for Z=l and multiplicity bin M3. 



The Streamer Chamber data fl^JIq] normally include all protons, whether free or bound in 
clusters. Generally the flow values using all isotopes oi Z = 1 and 2 are 10 — 20% larger 
than those for protons. In Fig. 3 the scaled flow values, Fs-, follow, within the uncertainties, 
a common trend with an initial steep rise and then an indication of a gradual decline. The 
Plastic Ball data are quoted with fairly small statistical uncertainties, ~ 4 — 7%, about twice 



the size of the data points in Fig. 3. For the Ar + Ph [0] and Ar + KCl Streamer 
Chamber data, we estimated the flow and statistical uncertainties from the vs. y data 
plots, for the appropriate multiplicity ranges. Our estimated uncertainties for these flow 
values are in the range of 15-23%. The 1800 A MeV Ar + KCl Streamer Chamber data, as 
analyzed in Ref. 7, produce a very large flow and scaled flow value [Fg — 93 MeV/c per 



nucleon which is off the plot scale and not included in Fig. 3). 

In summary, we have determined the nucleon flow for Ni induced collisions over an 
energy range of 400 to 1970A MeV. For these Ni systems flow, F, as measured by the 
change with rapidity of the average transverse momentum, rises with energy, and then 
declines. Comparison of our flow results with flow data from other mass systems is made 
by introducing a scaled flow, Fs = F/ {A\^^ + A^^^) which exhibits a nearly universal flow 
energy dependence. 
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FIG. 1. S curves for the four Ni + Au energies, using protons only in multiplicity bin M4. 

is the average of the x component of the momentum (see text), y' is the rapidity in the lab frame 
scaled by the rapidity of the beam. 
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FIG. 2. Proton flow as a function of beam energy per nucleon in multiplicity bin M4. The top 
panel is for Ni + Au EOS data. The solid symbols are the flow values from the graphs in Fig. 1 
while the open symbols are for soft and hard equation of state BUU Ni + Au calculations with 
an equivalent impact parameter distribution. The bottom panel consists of flow values determined 
from Ni + Cu EOS data. 
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FIG. 3. Scaled flow values vs. beam energy per nucleon for different projectile-target systems 

for Plastic Ball multiplicity bins 3+4. In the EOS and Plastic Ball data all isotopes of Z=l,2 are 

included. For the Streamer Chamber all free and bound protons are included. To improve the 

distinction between data points at the same beam energy, some of the beam energy values have 

been staggered around by as much as 20 MeV. 
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